1. Joint Modeling
1.1 Trace plots for convergence check
The current MCMC setting is:
- 100,000 iteration;
- 90,000 burn-in;
- 10 thinning.
1.2 Gelman and Rubin’s convergence check
## Potential scale reduction factors:
##
## Point est. Upper C.I.
## HDevsum 1.01 1.02
## LDevsum 1.02 1.05
## dh0 2.36 3.84
## dh1 5.65 11.83
## dh2 3.62 7.09
## dl0 1.62 2.49
## dl1 2.70 5.11
## dl2 2.89 5.27
## dl3 1.22 1.57
##
## Multivariate psrf
##
## 4.85
1.3 ACF Plots
Here we plotted ACF plots for the following variables:
- Total deviance;
- Variables that didn’t pass the convergence check.
1.4 WAIC results
| LevelH | LevelL | |
|---|---|---|
| DIC | 1184.41236 | 22680.2827 |
| DIC3 | 1139.45554 | 20662.2708 |
| PWAIC | 33.77812 | 948.6484 |
| WAIC | 1158.38855 | 21010.0069 |
2. Separate Modeling
2.1 Trace plots for convergence check
The current MCMC setting is:
- 100,000 iteration;
- 90,000 burn-in;
- 10 thinning.
2.2 Gelman and Rubin’s convergence check
## Potential scale reduction factors:
##
## Point est. Upper C.I.
## HDevsum 1.01 1.02
## LDevsum 1.02 1.06
## dh0 4.37 7.72
## dh1 7.31 12.81
## dh2 3.59 6.71
## dl0 1.42 2.11
## dl1 1.42 2.07
## dl2 1.12 1.33
## dl3 2.24 3.69
##
## Multivariate psrf
##
## 7.27
2.3 ACF Plots
Here we plotted ACF plots for the following variables:
- Total deviance;
- Variables that didn’t pass the convergence check.
2.4 WAIC results
| LevelH | LevelL | |
|---|---|---|
| DIC | 1377.23924 | 24118.439 |
| DIC3 | 1243.55301 | 21367.683 |
| PWAIC | 89.58375 | 1524.755 |
| WAIC | 1311.70933 | 22092.162 |